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ACRONYMS 

CAS Climate Adaptation Services 

CWM City water map 

DA Drought activity 

DEM Digital elevation model 

EAI Expected Annual Impact 

FLOPROS Flood Protection Standards 

GCM General circulation model 

GDP Gross domestic product 

GHG Greenhouse gas 

GTSM Global Tide and Surge Model 

IPCC Intergovernmental Panel on Climate 
Change 

IVM Institute for Environmental Studies 

MERIT Multi-Error-Removed Improved-Terrain 

PBL Netherlands Environmental Assessment 
Agency 

PDSI Palmer Drought Severity Index 

RCP Representative concentration pathway 

SPEI Standardized Precipitation and 
Evapotranspiration Index 

SPI Standardized Precipitation Index 

SRI Standardized Runoff Index 

SSMI Standardized Soil Moisture Index 

SSP Shared socio-economic pathway 

UNDRR United Nations Office for Disaster Risk 
Reduction 

WASP Weighted Anomaly of Standardized 
Precipitation 

US United States (of America) 
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Climate change, urbanisation and population growth are causing water availability and safety issues in cities around the world. Of all 
1 a natural hazards, flood and drought affected the most people between 2000 and 2019. Around 1.65 billion people were impacted by floods 
= and 1.43 billion by drought.’ Whether it be too much water too quickly, with extreme rainfall resulting in flash flooding in urban areas,? or 
Introd uction not enough water, because of long-term drought conditions;? no city is immune to these climatic changes. Sea-level rise and storm surges 
could force hundreds of millions of people in coastal cities from their homes at a total cost of more than US$ 1 trillion a year by 2050.4 
This increased migration and urbanisation is putting pressure on city water systems and infrastructure. Floods and droughts are, therefore, 

threatening the welfare of urban citizens and the critical services and operation of cities. 


Cities and their leaders need to understand this complex and imminent threat and its consequences. Knowing how much flooding and 
drought a city will experience and the potential impacts of these hazards will prompt city authorities and financiers to take action on 
adaptation and resilience. To this end, the analysis behind Water-safe Cities report quantifies the future hazards of flooding and drought in 
97 C40 cities by 2050 and the expected annual impact this will have on populations, urban infrastructure and the gross domestic product 
(GDP) of a city. The modelling and analysis have been used to develop interactive maps hosted on the Water-safe Cities webpage and the 
key headlines and messages of the research. 
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2. 


About this report 


This Technical Report details the methodology 
developed by the Institute of Environmental Studies 
and Climate Adaptation Services, in collaboration 
with C40 Cities, of the key communication tools 
behind the project. 


We first look to quantify the flood and drought 
hazard to C40 cities. Flood and drought risk are 
quantified in line with the common definition of risk 

as a function of hazard, exposure and vulnerability 
(see the Definitions section). Specifically, we assess 
three types of flooding (riverine, stormwater 
and coastal), as well as two types of drought 
(agricultural and hydrological). Because of data 
limitations, we have assessed stormwater flooding 
using the probability of extreme rainfall events as a 
proxy. To calculate the socio-economic impact for 
riverine and coastal flooding, we use state-of-the- 


PART I 


art global flood risk models, taking into account 
future socio-economic scenarios and using shared 
socio-economic pathways (SSPs) based on land- 
use projections from the 2UP exposure model.5 


We measure hydrological drought by estimating the 
annual deficits in urban water supply using modelled 
discharge time series and empirical extraction 
data. For agricultural drought, we assess modelled 
rootzone moisture shortfalls in agricultural areas 
around a city. After providing a socio-economic 
context for each city, we monetise the risks in 
US dollars for riverine and coastal flooding and 
hydrological drought using depth-damage curves 
and unit costs for water production/savings. 


The analysis is divided into two parts: 


A climate hazard analysis: populating our models with hydrological parameters to 
estimate the flood and drought hazard for 97 C40 cities. 


PART Il 


A socio-economic risk analysis: combining hazard, exposure and vulnerability 
parameters to calculate socio-economic risk and the probabilistic impacts of flooding 


and drought. 
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3. 
How we define flooding 
and drought 


Here, we define the climate hazards we used to produce our climate and socio-economic risk analyses and the interactive global maps. 


Flooding 
A flood is an overflow of water from a water 
conveyancing or storage system. Flooding is 
normally divided into three types: stormwater 
flooding (originating directly from rainfall), 
riverine flooding (originating directly from 
rivers) and coastal flooding (originating directly 
from the sea). Heavy rainfall that causes rivers 
to burst their banks is still considered riverine 
flooding, as the river is the source of the flood. 
We assess all three types of flooding in our 
analysis. 


Drought 

A drought is essentially a temporal shortage 
of water relative to a climatic baseline or 
historical average. However, water deficits 
can occur in different parts of the hydrological 
cycle, which makes a comprehensive definition 
challenging." In addition, droughts are often 
slow and creeping in nature, which makes it 
hard to define their onset, duration and end.? 
It is also hard to delineate the spatial border 
of a drought, as it is not directly observable.? 
Consequently, scientists and policymakers 
have put forward a variety of alternative 
definitions in recent decades, sparking the 
creation of more than 150 drought indices. 1° 


In recent decades, droughts have frequently 
been divided into three distinct types: 


1. Meteorological drought - a shortage of 
rainfall." 
2. Agricultural drought - a shortage of 


soil moisture that would otherwise be 
available for crop and vegetation 
growth. ? 

3. Hydrological drought - a lack of (sub-) 
surface water, including streams/rivers 
and roundwater.¥ 
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Thus, an appropriate definition of drought 
depends on the stage of the hydrological 
cycle. This, in turn, comes down to the 
(expected) impacts a drought will have on the 
area in question.^ In section 5.1, we assess the 
expected impacts of drought in urban areas 
and determine the type(s) of drought relevant 
to this project. 


A different but related concept is that of water 
scarcity, which is excess water demand relative 
to supply for a certain economic activity. 
Water scarcity differs from drought in that 
it is considered a systemic state (something 
permanent or ‘normal’) rather than a temporary 
phenomenon.” This study will focus primarily 
on drought, as it is more directly climate related 
than water scarcity, based on the definitions 
used here. For a practical example of water 
scarcity and its impact on cities, please see 
our accompanying web narratives. 

For detailed descriptions of how we define 
risk, please see the Definitions section. 
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4. 
Cities in this analysis 


Below are the 97 C40 cities that are part of this research. 
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Figure 1. C40 cities 
Abidjan Chengdu Guangzhou London Phoenix Stockholm 
Accra Chennai Hangzhou Los Angeles Portland Sydney 
Addis Ababa Chicago Hanoi Madrid Qingdao Tel Aviv 
Amman Ciudad de México Heidelberg Medellin Quezon City Tokyo 
Amsterdam Copenhagen Ho Chi Minh City Melbourne Quito Toronto 
Athens Curitiba Hong Kong Miami Rio de Janeiro Tshwane 
Auckland Dakar Houston Milan Rome Vancouver 
Austin Dalian Isanbul Montreal Rotterdam Venice 
Bangkok Dar es Salaam Jaipur Moscow Sao Paulo Warsaw 
Barcelona Delhi NCT Jakarta Mumbai Salvador Washington D.C. 
Beijing Dhaka Johannesburg Nairobi San Francisco Wuhan 
Bengaluru Dubai Karachi Nanjing Santiago Yokohama 
Berlin Durban Kolkata New Orleans Seattle Zhenjiang 
Bogota (eThekwini) Kuala Lumpur New York City Seoul 
Boston Ekurhuleni Lagos Oslo Shanghai 
Buenos Aires Freetown Lima Paris Shenzhen 
Cape Town Fuzhou Lisbon Philapelphia Singapore 
Guadalajara 
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5. 
Methodology 


Our overall methodology consists of the four steps 
outlined in Figure 2. 


1. Populating our models : 


Result: daily hydrological parameters : 


Result: Estimation of hydrological hazards for flooding : 
and drought 
: 


3. Calculating socio-economic impacts 
Result: socio-economic impacts for flooding and drought, 


for different return periods 


o 
Os 
e. 2. Developing global hazard maps === PARTI 
A ji 


^ PART II 
4. Estimating probabilistic risk : 

Result: average impacts per year (e.g. expected annual damage 

or affected people) for 97 C40 Cities 


+ Fluvial flooding City population 
* Coastal flooding Number of citizens | Energy 
> * Pluvial flooding i : Number of power generation facilities 
ve È Gross domestic product impacted by flooding 
©. ? | (GDP) 
PA Healthcare 
97 CITIES e eee Urban land-use Number of healthcare facilities 
ELE * Hydrological drought Residential, commercial, ' impacted by flooding 
and industrial 
2050 RCP45 & RCP85 
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5.1 
Part I: Climate hazard analysis 


For our climate hazard analysis, we review the first two steps in creating and populating global climate models to produce 


our global hazard maps. 


d* 


1. Populating our models 


Result: daily hydrological parameters 


2. Developing global hazard maps =- PARTI 


Result: Estimation of hydrological hazards for flooding 


and drought 


Figure 2a. Summary of the methodological approach (Part I: Climate hazard analysis) 


Any climate risk model requires comprehensive 
data and the integration of hazard, exposure and 
vulnerability parameters to achieve the best possible 
estimate of climate risk. In other words, the more 
detailed the input information on these factors, the 
closer the estimate will be to reality. Furthermore, 
a distinction is often made between models and 
studies that simulate only the climate risk of the 
current situation and those that estimate the climate 
risk of future scenarios. In our methodology, hazard, 
exposure and vulnerability can be statically (unable 
to evolve over time) or dynamically (able to evolve 
over time) added to the model. ' 


To make predictions, climate risk assessments often 
make use of representative concentration pathways 
(RCPs)” for scenarios on future greenhouse gas 
(GHG) concentrations and SSPs for scenarios on 
socio-economic development. RCPs and historical 
data on GHG concentrations are used in general 
circulation models (GCMs). A GCM models the 
climatic state of the earth based on a set of 
boundary conditions (mainly solar radiation and 
GHG concentration). The output of a GCM (for 
example, atmospheric temperature or rainfall) is also 
the input for hydrological models, which use this 
data to simulate hydrological processes, resulting in 
runoff/streamflow and groundwater recharge. 
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Next, the magnitudes and corresponding 
probabilities of possible hazardous events can be 
extracted. The probability of a hazardous event 
happening is often expressed in conjunction 
with a 'return period', especially in the case of 
flooding. This is the average amount of time 
(usually in years) between two such events. It 
is more precisely defined as the inverse of the 
probability that an event is exceeded in any one 
year: Rp = 1/P. 
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Understanding return periods 

Return periods are often determined using 
extreme value analysis. This involves conducting 
an assessment based on past events above a 
certain threshold or the maxima/minima each 
year. This series of events is ranked from the 
most extreme to the most common, then a 
statistical (extreme value) distribution is fitted to 
the data. From this fit, it is possible to derive the 
expected magnitude of an event for each return 
period. It is also possible to extrapolate this when 
acknowledging the uncertainty of exceeding 
known events. Hence, a 10,000-year event can 
be estimated without ever being observed, albeit 
with a high uncertainty interval. 


Rp is the return period in years and P is the 
probability of an event being exceeded in any one 
year (exceedance probability). For example, if a 
flood of volume X happens statistically once every 
100 years, the return period of that flood is said to 
be 100 years. The probability of a flood of at least 
that magnitude occurring in a given year becomes 
(1/100)*100% = 1% a year. Obviously, these return 
periods tend to be based on historical data; the 
exact statistical probability of any hazardous 
event occurring is unknown. This also means that 
a 100-year event might not occur in the next 100 
years or that a 10O-year event might occur several 
times in 100 years. For instance, high water levels 
and consequent flooding happened along the 
Rhine and Meuse rivers in both 1993 and 1995. 
Statistically, these events were estimated to occur 
only once every 25-50 years. '? 


Note: We use return periods in this study to 
indicate the probability of flooding, but do not 
use them in our drought risk assessment. 


It should be stressed that return periods are 
dynamic, so evolve over time. A return period 
is linked to a certain magnitude of event. If, due 
to climate change, events of this magnitude 
happen more frequently, their return period will 
decline. For example, a 1O-year return period 
rainfall event might equal a downpour of 40mm/ 
hour in a certain city. Hence, every year, there is 
a (1/10)*100% = 10% chance that this will happen. 
Should precipitation patterns change due to 
climate change and that 40mm/hour event now 
have a 20?6 probability of occurring per year, this 


would alter the return period to 1/(20%/100%) 
- 5 years. Thus, the return period for the event 
would halve because the frequency of the event 
had doubled due to climate change. 


Another general assumption is that the Global 
Human Settlement Model Grid (GHS-SMOD) 
city border dataset is the most suitable dataset 
for global flood and drought assessments. 
Countries define their cities in different 
ways,? which makes it difficult to consistently 
compare urban flood and drought risk between 
countries. This can be overcome by choosing 
a physical definition for cities rather than an 
administrative one? The GHS-SMOD applies 
such a definition and bases it on building 
and population density thresholds,” making 
it consistent between countries.? However, a 
physical definition also means that adjacent 
cities can be merged into one agglomeration 
within the dataset. Still, it is likely that impacts 
in one city will spill into adjacent cities, as they 
are probably economically connected. 


WATER SAFE CITIES | TECHNICAL METHODOLOGY 


Step 1: 
Populating our models 


For this study, we use data from the PCRaster 
GLOBal Water Balance model (PCR-GLOBWB)?5 
as the basis for simulating river flooding and 
the water supply side for hydrological drought. 
In this study, we used the new version of 
the model from 2018.24 PCR-GLOBWB is a 
global hydrological model that simulates daily 
hydrological parameters such as discharge and 
the water volume that resides in the floodplains 
of rivers at a resolution of 5’ x 5’ (around 10 km 
x 10 km at the equator). 


The output of the climate model is used directly 
or indirectly to calculate hazard indicators. For 
our purposes, indirectly means that a global 
hydrological model is used to calculate the 
corresponding values of discharge and rootzone 
moisture (the volume of moisture in the top 1 
metre of the soil) from these climate model 
outputs. These values are then inserted into 
another model that calculates the hazard. 


Riverine flooding 

For riverine flooding, PCR-GLOBWB overbank 
flow volumes are used to develop inundation 
maps using the mass-conservative flooding 
algorithm developed by Winsemius et al. (2013). 
We use this to generate global maps showing 
the extent and depth of flooding for different 
flood return periods (2-, 5-, 10-, 25-, 50-, 100-, 
250-, 500- and 1,000-year return periods) at a 
horizontal resolution of 30” x 30” (around 1 km 
X 1 km at the equator). 


Coastal flooding 

We model coastal flooding by (a) modelling 
the sea level for different return periods and 
(b) simulating the on-land inundation from 
these sea levels. The sea level comprises three 
components: the tide, the storm surge and wind 
waves. We do not include the latter here. For 
step (a), we use the Global Tide and Surge Model 
(GTSM) to simulate the sea levels resulting from 
storm surge and tide for the return periods 
as described for riverine flooding. We use the 
latest version of GTSM.? The model resolution 
is 2.5 km from the coast (1.25 km in Europe). 
GTSM is forced with ERAS high-resolution 
meteorological forcing data. This dataset has 
a spatial resolution of 30 km x 30 km, enabling 
it to better resolve tropical cyclones than earlier 
versions of GTSM$ The baseline period for these 
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simulations is 1979-2017 (for climate forcing), 
with a mean sea level equivalent to the 1986- 
2005 period. For future simulations, probabilistic 
sea-level rise projections were carried out using 
the method presented by Le Bars et al. (2017), 
based on CMIP5 data for RCP 4.5 and RCP 8.5. 


Next, in step (b), the extreme sea levels for 
different return periods are converted into flood 
hazard maps using the approach described in 
Tiggeloven et al. (2020). This approach uses a 
GIS-based routine, in which the GTSM extreme 
sea levels are projected to the nearest coastline, 
with inundation taking place in areas that are 
hydrologically connected to the sea at that 
extreme sea level. The model uses the Multi- 
Error-Removed Improved-Terrain (MERIT) Digital 
Elevation Model as underlying topography (30” 
x 30”) and we use a resistance factor to simulate 
the inward reduction of flooding land, as tides 
and storm surges have a limited timespan. 


Stormwater flooding 

For stormwater flooding, to represent the future 
hazard, we use daily rainfall data that come 
directly from GCMs.?6 We used rainfall data from 
ISIMIP2b for the following GCMs: GFDL-ESM2M, 
HadGEM2-ES, IPSL-CM5A-LR, and MIROCS5, 
again for RCPs 4.5 and 8.5. 


DROUGHT 


Because of their complexity and diversity, 
it is often difficult to find a comprehensive 
set of indicators for drought risk. This has led 
to a variety of approaches and datasets in 
global drought risk assessment. Most drought 
indices focus on rural rather than urban areas, 
and drought risks to or impacts on cities are 
investigated more often on a local or regional 
scale.” Moreover, drought risk damages have 
not been widely monetised on a global scale 
before.?? 


As mentioned, there is no single, comprehensive 
drought indicator. Whether a given indicator 
is relevant to this assessment depends on the 
water sources affected by drought and whether 
this has an impact on cities. We, therefore, 
gathered information on drought impacts in 
urban areas 


to determine which parts of the hydrological cycle 
were relevant. We then chose two drought indicators 
for this assessment: hydrological drought and 
agricultural drought. We elaborate on this process 
of indicator selection below, based on their impacts 
on urban areas. 
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Figure 3. 


Left: Water resource dependency per C40 city, 
according to the City Water Map database (CWM)? 


Urban water sources 

Drought impacts in cities are most often related to 
urban water supply and agricultural practices. To 
assess shortfalls in urban water supply (domestic 
and industrial), it is important to identify the source 
of this water. McDonald et al. (2014; 2016) created 
a dataset in which they showed the municipal 
water sources for 534 large urban agglomerations 
(with a population of more than 750,000) globally. 
This dataset also includes the source type (for 
example, surface water or groundwater), as well as 
an estimate of the volume of water withdrawals for 
several agglomerations. According to this dataset, 
cities most often depend on either groundwater or 
surface water. A city is deemed reliant on a specific 
source if 50% or more of its water originates from 
that source. Sporadically, other water sources (such 
as desalination) are used, but these do not meet the 
50%-plus criterion for any city. Moreover, saltwater 
bodies are generally too large to be affected by 
drought and hinder desalination. Hence, a drought 
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in surface water or groundwater is likely to have the 
greatest impact on an urban area in terms of public 
water supply. This is what is called a hydrological 
drought. The water dependency of the C40 cities 
can be found in Figure 3 (spatially) and Table 1 
(overview). 


Volume of water withdrawal per C40 city 


ra ae í 4, 


Sources & ownership: 

Data frome C40-Gities; City Water Map VLG (McDonald et al., 2016); Natural Earth Data (background map). 
Created with QGIS 3.10.8 

Map created by: Tristian Stolte (IVM VU Amsterdam) 


Right: Volume of water withdrawals per C40 city, 
according to the CWM 


SELECTING INDICATORS 

Establishing two main categories of drought impact 
in urban areas calls for two different drought indices, 
each able to estimate the drought hazard for one of 
those categories. 


Spatially standardised indices are often preferred 
for studies on a global scale because: 


e They allow a comparison of values between 
different climatic zones. 

* They require fewer data points. 

e Their interpretation is intuitive. 

e They perform robustly over multiple studies. 


This type of indicator is especially useful when 
no absolute threshold is known beforehand. The 
threshold below which drought occurs will then 
be based on local climatic conditions and the 
corresponding statistical occurrence of dryer-than- 
normal events. 
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Percentage of C40 
cities 


Table 1. Water-source dependency of C40 cities (surface 
water and groundwater) 


Dependent on 


Note: Not all 97 cities are included in the water resources 
database and there are cities without estimates for 
water withdrawal volumes. An alternative methodology 
for classifying hydrological drought from a qualitative 
methodology methodology in Annex 2. 


HYDROLOGICAL DROUGHT 

We came up with a non-standardised indicator for 
the urban water supply in this study as it enabled 
an absolute threshold to be derived for each water 
source. This indicator compares water availability 
(water supply) from the hydrological model 
with estimates of daily water withdrawals (water 
demand) per water utility. The latter were found in 
the CWM database,’ which contains the locations 
of water utilities for at least 88 of the 97 C40 cities. 
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However, there are only withdrawal estimates for 76 
of the cities. estimates of daily water withdrawals 
(water demand) per water utility. The latter were 
found in the CWM database,’ which contains the 
locations of water utilities for at least 88 of the 97 
C40 cities. However, there are only withdrawal 


estimates for 76 of the cities. Where we were not 
able to account for the remaining C40 cities we 
have provided a qualitative methodology in Annex 
2. 


Furthermore, drought volumes are standardised 
per 1,000 inhabitants. Population counts are 
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also extracted from McDonald et al. (2016). Thus, 
the hydrological drought hazard indicator is the 
expected annual volume of drought per 1,000 
inhabitants. 


AGRICULTURAL DROUGHT 

The second main impact category is agriculture, 
which is caracterized by two drought types, 
depending on which irrigation strategy is chosen: 
either soil moisture drought for rainfed agriculture 
or hydrological drought for artificially irrigated 
agriculture. 


Agricultural drought has been selected due to 
its indirect impact on urban areas. Furthermore, 
agricultural drought has been used to supplement 
our hydrological drought analysis as it was not 
possible to undertake a complete analysis for 
several C40 cities. 


We use the Standardized Soil Moisture Index (SSMI) 
to assess agricultural drought, as it relates to rainfed 
agriculture, the most common type of agriculture.*5 
The SSMI calculates the statistical distribution of 
soil moisture volumes for each month in a given 
year based on time series of soil moisture for that 
month. A threshold of -1 standard deviation is used, 
which means that if soil moisture drops below one 
standard deviation of the mean value, we consider it 
a drought. is the expected annual volume of drought 
per unit of area. 


WATER SAFE CITIES | TECHNICAL METHODOLOGY 


Step 2: 
Developing global 
hazard maps 


We use the aforementioned hydrological and 
climate parameters to develop the indicators of 
flood and drought hazards mentioned. These, in 
turn, provide hazard information for the 97 C40 
cities. This section describes the development of 
each hazard map individually. Table 2 summarises 
the hazards and indicators used. To understand the 
resolution of each hazard map, please refer to the 
Resolution Tables in Annex 4. 


Flooding 


On a global scale, the flood hazard for a certain 
geographical area is generally represented by the 
inundation depth and extent corresponding to a 
specific return period (see Understanding return 
periods on page 9). In the case of riverine flooding, 
the inundation depth and extent are derived from 
the spatial propagation of water from the overbank 
flow corresponding to the relevant return period. 
For each pixel on the gridded inundation maps, we 
multiplied the inundation depth by the area of that 
pixel (the extent) to get the inundation volume per 
pixel. The total inundation depth is the sum of all 
pixel volumes inside the city. We then standardise 
this by dividing by the total city area to obtain the 
hazard indicator for riverine flooding in terms of 
volume per area. We use the same process to derive 
the coastal flooding hazard, except that the flood 
water comes from the overtopping of coastlines 
rather than riverbanks. 


The same approach cannot be used for stormwater 
flooding, as the model is not yet able to create global 
inundation maps for rainfall. Therefore, the change 
in extreme precipitation events per pixel is used 
as a proxy for future stormwater flooding hazard. 
The precise indicator is the factor change in the 
frequency of occurrence of a current precipitation 
event with a 10-year return period. For example, if 
a 10-year return period rainfall event of 40mm per 
hour becomes a five-year return period event in 
future, the hazard indicator takes on a value of 2, as 
the event becomes twice as frequent. 


Riverine and coastal flooding 

For riverine flooding, we produced hazard maps of 
the flood volume (the extent of the flood multiplied 
by the depth of inundation) for each city for different 
return periods, as described in section 5.3. 
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We took an ensemble mean of flood volumes per 
pixel, or the average flood volume for all GCMs. 
Next, we aggregated the flood volumes per city 
and standardised them to m3/km* for comparable 
results between cities. This gave us an estimate of 
the riverine flooding hazard per city for each return 
period and climate scenario (that is, historical/ 
current, RCP 4.5 and RCP 8.5). 


We derived the coastal flooding hazard map in a 
similar way, giving us an estimate of the coastal 
flood hazard in m?/km* for each combination of 
city, return period and climate scenario. The only 
difference is that coastal flood inundation maps are 
based on probabilistic sea-level rise projections, so 
are not calculated for different GCMs, but already 
statistically aggregated over different climate 
models. 


For the surface area, we delineated each city using 
the city borders from the GHS-SMOD.** 


In this study, we used a state-of-the-art method 
for our riverine and coastal flood assessments. 
The main limitations stem from the trade-offs 
in data availability and accuracy needed for 
global model analyses. We explain them further 
here, based on the work of Ward et al. (2020b) 
and Tiggeloven et al. (2020). The flood hazard 
modelling approach uses generic return-period 
maps, so does not consider specific flood 
events. This means, for instance, that a 1- 
in-100 event could correspond to different 
flood volumes in different parts of a city - and 
this is likely to be the case in reality. 
Moreover, the flood inundation simulations do 
not use a hydrodynamic modelling scheme 
and are carried out at a spatial resolution of 
30” x 30”, which is not sufficiently granular to 
include local topographical features in cities. 


Another limitation is that coastal and riverine 
flooding are considered independently of each 
other. This means, for instance, that the flood 
hazard simulations do not include downstream 
backwater effects in coastal cities. We, therefore, 
assume that rivers in coastal regions can flow 
away into the ocean without being impeded by 
high coastal water levels in deltas and estuaries 
(for example, caused by storm surge and high 
tide). The exposure data have a resolution of 
30” x 30” and do not differentiate between 
different kinds of urban fabric. However, we 
have made an estimate of the percentage split 
of urban area between residential and 
commercial. 


Also, standard depth- 
damage functions are used around the globe due 
to a lack of spatially disaggregated information. 


These curves only apply to direct damages and 
exclude indirect damages. 


Stormwater flooding 

For stormwater flooding, we calculated changes 
in the occurrence of extreme rainfall events. We 
chose the current 10-year return period event as our 
baseline and calculated how much the frequency 
of this event might change in future scenarios. As 
mentioned, return periods are dynamic. For each 
pixel inside a city, we took the total time series of 
daily rainfall from the historical data (1971-2006). 
For each year, we looked for the maximum rainfall 
volume per day and then sorted those maximum 
levels from large to small. We then searched for the 
10-year return period rainfall event and identified 
the corresponding rainfall volume. 


We performed the same sorting process on the 
future scenario data (2030-70) and identified the 
rainfall volume of the historical/current 10-year 
return period event in the future time series. We 
then calculated the corresponding future return 
period and divided it by the historical return period 
to obtain the change in frequency of occurrence. 
We then averaged these numbers over all climate 
models and all pixels inside a city to get an estimate 
of the stormwater flooding hazard per city per RCP. 
Cities were again delineated using GHS-SMOD. 


As mentioned, there is a big data challenge in 
assessing stormwater flooding on a global scale. 
This means we could only provide a simple proxy for 
the stormwater flood hazard here. Consequently, we 
have not estimated the damage from stormwater 
flooding. Nonetheless, the hazard maps highlight 
several cities as potential hotspots of stormwater 
flood risk, and this could be an incentive for more 
detailed studies on those cities. 


Furthermore, we have analysed the change in 
frequency of a 10-year rainfall event rather than the 
change in volume, as the non-granular resolution of 
the input data did not allow us to accurately estimate 
total rainfall volumes at city level. Hence, the change 
in frequency is probably more consistent over 
different spatial scales than the change in volume. 
The 10-year event itself was chosen based on the 
temporal resolution (daily) and temporal extent 
(35-40 years) of the input data. Higher-intensity 
events require longer time series, as they occur less 
often than lower-intensity events. 


WATER SAFE CITIES | TECHNICAL METHODOLOGY 


21 


Global gridded dataset 
of climate simulations for 

; current and future climate 
produced using the mean 
of five GCMs 


Global gridded dataset 
of climate simulations 


for 


current and future 


climate (GCMs) 


" 


P 2 


methodology behind riverine flooding 


HYDROLOGICAL MODEL Estimate overbank flow for nine 
different return periods 


PCR-GLOBWB mm » (ranging from a 
Hydrological model one-in-two-year event to a 
—— > estimates river discharge for one-in-1,000-year event 


current and future climate 


f it Y Fluvial flooding calculated for 
[XM \ 97 C40 Cities 


Inundated volume (m?/km7?) 


CLIMATE SCENARIOS 
2050 RCP4.5 & RCP8.5 
Estimate inundation volumes for 
nine different return periods for [RETURN PERIOD | 
all 97 cities chosen in the maps 
& 1in a 100 year 
a A. 97 C40 city borders INS 
E based on the GHSL — > 
< Depth 
Extent 


Coastal flooding 


Inundated volume (m?/km7?) 
for different return periods, 
for each city 


Digital Elevation Model 


MERIT 
Global Tide and Surge Model (GTSM) i 
simulations of sea levels resulting Estimate overtopping of £, 
from storm surge and tide for current coastlines for nine different 2050 RCP4.5 & RCP8.5 
and future climates return periods (ranging from a 
one-in-two-year event to a 


ic NE one-in-1,000-year event) chosen in the maps 
lin a 100 year 


Estimate on-land inundation volumes for 


nine different return periods for oo 
all 97 cities SO 
€ ae 
» | 97 C40 city borders ja oo 
^ | based on the GHSL yO” 
be Depth x 
300 
. : Extent tne metn 
ethodolo9y behind coastal flooding «mat of 
em cual su 
ay 9 = Figure 4€ visu 
exe" i \ 
qu? B For each pixel 
xe bo: Global gridded dataset determine one-in-10-year downpour 
A) of climate simulations "T (which isa rainfall event with a 
tor current and future Select total precipitation probability of 10% each year) 
climate (GCMs) per day for the global grid 
l Stormwater flooding 
"h increase in probability of a 
" current one-in-10-year 
B Analyse the change in frequency i 
estimate how often this downpour downpour for each city 
might occur in 2050 
dih iih 2050 RCP4.5 & RCP8.5 
imi iam d 
& m => — — um 
€ — 97 C40 city borders Eu X X he chosen in the maps 
a based on the GHSL = x = lin a 10 year downpour 
NOW 2050 
22 WATER SAFE CITIES | TECHNICAL METHODOLOGY 


Hydrological drought 

For hydrological drought, we created a time series 
of discharge per calendar month and calculated the 
average discharge per pixel for each month in the 
series. This resulted in 12 average discharge values 
per pixel, one for each month of the year. We then 
analysed the complete discharge time series and 
identified when the discharge was below average 
for the corresponding month. We assumed that 
shortages occurred below this threshold. We then 
added up the annual water shortages per city 
surface water source and per city. Note that several 
cities extract their water from multiple surface 
water sources. Next, we averaged the total water 
shortages per city per year over the climate models 
and divided this by the number of city inhabitants. 
This hazard is thus expressed per 1,000 inhabitants 
per city per RCP. Note, too, that this metric inherently 
includes exposure by considering the number of 
citizens. We used only the CWM population counts 
for the 2010 base year to separate out the climate 
effects on future changes in drought hazard. 


Lower range| Upper range 
Low 
Medium-low 


Medium 


Medium-high 


Table 2. Hydrological drought risk analysis for CAO 
cities 
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In summary, the main assumptions used in 
our hydrological drought hazard calculations 
likely lead to an overestimate of actual drought 
volumes from surface water. These could be 
either alleviated or aggravated by drought 
volumes that originate from groundwater 
sources. Moreover, these drought volumes 
potentially hit poor citizens harder than 
wealthier ones; we have not been able to 
account for this in our analysis. 


In addition, we also assume that cities use the 
water sources set out in the CWM database.“ 
Unfortunately, not all C40 cities are included in 
the database and the information is incomplete 
for some cities. For example, Cape Town uses 
several surface-water sources for its water 
supply?9 but these are not present in the 
CWM. Still, the database is the most complete 
overview of water withdrawals available on 
a global scale. For missing cities, we have 
provided a qualitative analysis and ranked their 
risk of hydrological drought below. A detailed 
summary of how this qualitative analysis was 
conducted can be found in Annex 3. 


Cities 
None 


Zhenjiang, Hangzhou, Shenzhen 


Berlin, Fuzhou, Rome, Madrid, Dubai, Kuala Lumpur, 
Milan, Singapore, Medellin, Stockholm, Abidjan, Quezon 
City, Ekurhuleni, Jaipur, Miami 


Lima, Accra, Beijing, Curitiba, Delhi, Jakarta, Nairobi, New 


York, Shanghai Washington, D.C. Johannesburg, 


London, Sydney, Cape Town, Melbourne, Rio, 
Copenhagen, Hanoi, Durban (eThekwini), Seoul 


Los Angeles, Santiago, New Orleans 


Agricultural drought 

For agricultural drought, we first used the SSMI to 
analyse time series of historical rootzone moisture 
and derive a standard normal distribution of rootzone 
moisture values per calendar month per grid cell. 
We identified a -1 standard deviation of historical 
rootzone moisture volumes for each grid cell, which 
we then translated back into an absolute volume of 
water. We assumed that there was a crop-damaging 
drought event if this threshold was exceeded, both 
in historical and future time series. We added up the 
total annual volume of water below the threshold for 
each grid cell. We then aggregated this information 
into one or more water provinces for each city, 
assuming that each city obtained its food from the 
water province(s) in which it was located. We also 
used physical agricultural area estimates from the 
MapSpam project to correct for the percentage of 
agricultural area in each water province. The last 
step in deriving the drought hazard for agricultural 
areas was to divide by the total surface area of the 
intersecting water provinces for each city. As we 
took the ensemble mean once again, the hazard 
was thus defined as the agricultural drought per city 
per climate scenario. As with urban water supply, 
we included exposure by defining the areas where 
potential damage could occur, namely, agricultural 
fields andthe cities to which they belonged. However, 
we did not directly include the exposed elements 
in the cities themselves (for example, population or 
companies affected), as this information could not 
be ascertained from the data in hand. 


In general, drought is very difficult to frame, 
especially on a global scale, so assumptions are 
required to quantify drought risk. In calculating 
occurrences of drought, a -1 standard deviation 
of historical rootzone moisture volume for each 
grid cell is assumed to be a crop-damaging 
drought event if this threshold is exceeded in 
either historical or future time series. 


Another key assumption is that crop failure 
in a water province surrounding a city has 
an impact on the city’s food supply and 
corresponding food prices. The import and 
export of food is not considered here due to 
data constraints, but it is likely that a city is not 
fully dependent on the farmlands in its own 
water province. However, even if crop failure 
in a region does not lead to food shortages in 
nearby cities, that region is still likely to suffer 
from local economic problems that may have 
chain effects for cities (for example, through 


24 


food distribution and processing). Hence, 
this assumption does not, in and of itself, 
overestimate the drought costs for a city. 


If a water province is shared by two or more 
cities, the drought hazard from that water 
province is shared equally between them. In 
reality, some cities may depend more heavily 
than others on the crops in a water province, 
so will experience a greater share of the total 
drought risk. Similarly, all crops are assumed to 
be equally affected by drought, although some 
may be more resistant than others. This means 
that some cities that depend on drought- 
resistant crops will experience a lesser drought 
impact than others. Moreover, the spatial 
resolution of 350 x 30 is not sufficiently 
granular to resolve agricultural droughts in 
very small water provinces. 


On a technical note, the SSMI-based drought 
threshold is only calculated for those pixels 
that have at least 8096 of the timesteps in the 
data. This cut-off point was chosen to exclude 
thresholds defined on just a few data points. 
The 8096 is a somewhat arbitrary number 
based on the authors' expert judgement. 
Furthermore, the SSMI threshold was set at -1, 
which detects a moderate drought or worse 
according to Komesku (1999) or a moderate to 
severe drought according to the classification 
of Agnew (2000). The one-month lag time was 
chosen based on the assumption that cities 
needed to deal with drought immediately 
when it occurred and that no water was stored 
or saved in between. Moreover, if we used a 
larger time lag when calculating the annual 
deficit, this would result in double counting. 
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Summarising our climate hazard analysis 


The table summarises the climate hazards analysed in this research and the indicators used for each hazard, 
as well as the unit in which we produced results, as displayed in the Water-safe Cities interactive maps. 


Climate model 


Indicator 
ensemble 


Spatial scale 


Riverine flooding Total flood volume Yes, averaged Per city 


Yes, inherently 


Total flood volume included 


Coastal flooding Per city 


Change in the frequency of a 
Stormwater flooding one-in-10-year precipitation 
event 


factor of change; 


Yes, averaged 5 i 
dimensionless 


Per city 


Per city 


Hydrological drought Urban water supply shortage Yes, averaged m3/year/1000 inhabitants 


(per inhabitant) 


Agricultural drought SSMI Yes, averaged 


Per water province m3/kmz2/year 


Table 3. Summary table of the climate hazards analysed and indicators used 


27 


5.2. 


Part Il: Socio-economic risk analysis 


The global hazard modelling methodology described in section 5.1., Part I: Climate hazard analysis provides 
the basis for Part IIl: Socio-economic risk analyses, in which we model the impacts and risks. We explain the 


methodology behind Part II later in this chapter. 


different return periods 


affected) for 97 C40 Cities 


4. Estimating probabilistic risk 


Result: Average impact per year (for example, expected annual damage and people 


3. Calculating socio-economic impacts 


Result: Socio-economic impacts for flooding and drought for 


-— PARTII 


Figure 2b. Summary of the methodological approach (Part Il: Socio-economic risk analysis) 


Indicators of exposure and vulnerability 


To understand the impact risks to cities from flooding and drought, we first state the indicators used to factor 
in exposure and vulnerability. For a full definition of how we calculate risk, please see the Definitions. 


This analysis focuses on the global scale, requiring 
detail to be weighed against the possibility of 
gaining a larger, consistent overview of flood and 
drought risk for the 97 C40 cities. Our approach 
makes use of state-of-the-art global-scale models, 
resulting in first-cut cost estimates. These outline 
the potential impact a drought or flood could have 
on each city, so can be used to raise awareness. It 
should be noted that our estimates do not include 
all potential flood and drought impacts. For 
example, we only examine direct impacts and do 
not look at the indirect impacts or ‘chain effects’ 
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that these hazards can cause. Such chain effects 
are often event and location specific, and methods 
to assess them are still difficult to implement in 
global disaster risk models, as they need highly 
accurate input data.*?? Moreover, we do not examine 
the potential impacts of floods and drought across 
all potentially impacted sectors, as we would have 
to take into account all chain effects to do so. 
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Flooding 


Flood exposure centres on three metrics: (1) 
population, (2) GDP and (3) urban land use. We 
estimate the first two using global datasets for 
population and GDP, each with a spatial resolution 
of 30” x 30”.4° These datasets include forecasts of 
the current situation and future socio-economic 
scenarios (SSPs) in the form of total population 
estimates per grid cell and total GDP per grid cell 
in 2005 US dollar terms.^ We chose this price 
unit because it is commonly used in disaster risk 
studies,? thus facilitates cross-study comparability. 
We modelled the urban land-use indictor using 
the 2UP model developed by the Netherlands 
Environmental Assessment Agency (PBL).^ This 
model is essentially an urban/non-urban land-use 
dataset that estimates the percentage of build-up 
in urban areas. Every urban land-use cell in 2UP has 
a fixed mix of three urban categories - residential, 
commercial and industrial - and each class is 
assigned its own maximum potential damage level. 
The spatial resolution of the data is 30” x 30” and 
the data are available for the current situation and 
future scenarios using SSPs. 


We express flood vulnerability using depth-damage 
curves, which depict the relationship between a 
certain depth of inundation (water level) and the 
expected costs for a certain type of building and/or 
location. In essence, a flood with a low inundation 
depth causes less damage than a flood with high 
inundation depth.^^ This is predominantly associated 
with the urban land-use exposure metric. 


Infrastructure exposure 

We also use two other metrics to denote the global 
(change in) vulnerability under different climate 
scenarios: (1) the number of power stations affected 
and (2) the number of hospitals and healthcare 
facilities affected. These are both defined as the 
total number of power stations/hospitals located 
in inundated grid cells. Taking global datasets for 
hospitals*® and power stations and overlaying them 
with our inundation maps, we selected the hospitals 
and power stations located in a given city. We then 


We did not consider local flood adaptation 
measures (such as elevated building sites 
or flood walls), so not all exposed stations 
or hospitals are likely to be affected. These 
buildings are units of critical infrastructure, 
however, and flooding can have relatively far- 


reaching effects on them and the cities they 
serve. Therefore, these metrics are intended 
to show the potential global number of critical 
infrastructure units within the exposed area as a 
measure of the overall vulnerability of the cities. 


overlaid the hospitals and power stations with the 
flood hazard maps. We then counted the number 
of hospitals and power stations^? in the inundated 
zone (inundated zone - all pixels that had a flood 
volume » O) and noted their ID numbers/names. 
We did this for every flood hazard map, so for 
every combination of climate model (GCM), climate 
scenario (historical, RCP 4.5 and RCP 8.5) and return 
period. 


Conventional exposure metrics used in flooding 
do not fit with the current approach to calculating 
drought risk in cities. The elements exposed to 
drought are inherently included in the hazard 
analysis for this study. For agricultural drought, 
the exposed elements are the agricultural areas in 
the water province(s) surrounding a city, while for 
hydrological drought, the exposed elements are the 
surface water-source locations and city populations. 
In this sense it is not possible to state the number 
of people in a city who are exposed to both 
hydrological and agricultural drought, as the city 
as a whole experiences drought, so the population 
within the city does as well. 


Drought vulnerability is not inherent in the hazard 
analysis, but included qualitatively after the socio- 
economic impact estimates in the Water-safe 
narratives. Vulnerability is represented by several 
metrics, depending on data availability and the link 
between the drought hazard and its impact on the 
exposed elements. These indicators are listed in 
Table 4. 


30 


p unn Reese m 
mme LE 


WATER SAFE CITIES | TECHNICAL METHODOLOGY 
^ nig: Sonic m RR 


Pm ~~ 


Step 3: 


Calculating socio-economic impacts 


Flooding 


The ingredients of the hazard, exposure and 
vulnerability metrics described in the previous 
sections could now be integrated into a socio- 
economic impact risk) estimate for each of the 
C40 cities. We overlaid the flood hazard maps 
for the different combinations of GCM, RCP and 
return period with population and GDP maps for 
the different SSPs. This resulted in a set of maps 
encompassing th e affected population and GDP 
per grid cell for each unique combination of GCM, 
RCP, SSP and return period. We also overlaid the 
urban land-use map for the different SSPs with 
the flood hazard maps. The percentage of urban 
built-up was then subdivided into commercial, 
residential and industrial build-up areas and linked 
to their corresponding depth-damage curves. 
These damages were summed to get total damage 
per grid cell. This gave us urban damage per GCM, 
RCP, SSP and return period combination. We then 
added up the impacts for each combination of 
GCM, RCP, SSP, return period and impact type for 
each individual city. Note that GCMs are relevant to 
riverine flooding but not coastal flooding. 


We estimated drought costs by calculating the cost 
of freshwater production/saving as a proxy. The 
need to use a proxy stemmed from the absence of 
an applicable method for calculating drought costs 
in urban regions on a global scale. Drought can have 
myriad (in)direct impacts that are hard to quantify, 
as described previously. The rationale for using 
the cost of freshwater production/saving was that 
this would quantify the cost of solving a drought 
of volume X, by producing X amount of water to 
maintain business as usual and, thus, prevent any 


Adaptation measure 


Increase reservoir storage capacity 


Reuse urban residential/industrial water Variable between water provinces 


Variable between water provinces 


other impacts from occurring. Although these costs 
are based on adaptation measures, we do not regard 
them as the true costs of adaptation. The measures 
chosen are all primarily mitigating solutions: they 
do not prevent drought, nor are they long-term, 
sustainable solutions. 


Freshwater production cost/savings are based on 
three adaptation measures for which global unit 
cost estimates are available from the water scarcity 
literature." These measures are: (1) the increase in 
reservoir storage capacity, (2) the reuse of urban 
residential/industrial water and (3) the increase in 
desalination (Table 4). 


Increased reservoir capacity costs are estimated 
based on Ward et al. (2010a). They vary from 
water province to province and are based on the 
relationship between mean basin slope and the 
unit cost per mž for 11 different reservoir sizes. Each 
water province used in this study was assigned a 
value from a water province in Ward et al. (2010a), 
based on the largest overlap, because the shapes 
occasionally differed from the more recent version 
of the report used here Straatsma et al. (2020). 


There are two elements to the costs of desalination: 
a constant value of US$ 1/m? conversion cost and a 
variable price for water transport between the city 
and the nearest coastline.^? Coastlines are obtained 
from Natural Earth and come from the coastline 
dataset V4.1.0.°° The transportation costs are divided 
into a horizontal and a vertical component, which 
carry a unit cost of US$ 0.05/m* per 100 km and 
US$ 0.05/m? per 0.1 km, respectively (see Table 4). 
Horizontal transportation distances are determined 


Unit cost (US$ 2005/msz) 


Ward et al. (2010a) 


Straatsma et al. (2020) 


1 (constant component) + 


Increase desalination 


0.05 per 100 km horizontal or 
0.1 km vertical transport (variable 


Straatsma et al. (2020); Zhou & 
Tol (2005) 


Table 4. Global average unit cost of freshwater production as found in the cited literature 


31 


DROUGHT 
~ 


pd 


Water deficit 
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NO 
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cheapest option 


YES 


Solve 5096 of the drought with the 
reuse of urban residential/industrial 
water and check for the next cheapest 


option 


Solve the remainder of the 
drought with the next cheapest 
option and aggregate the costs 


Figure 7. Process for calculating the cost of hydrological drought 


by measuring the shortest linear distance 
between a city's boundaries (as derived from 
GHS-SMOD) and the nearest coastline, using 
QGIS 3.10.11. The vertical transportation distance 
is the sum of all absolute height differences 
between consecutive pixels from the Shuttle 
Radar Topography Mission Digital Terrain Model 
(SRTM DTM)?! V4.1 that intersect with the 
horizontal transport line. The SRTM DTM has a 
resolution of 3 arcseconds (3). Lastly, we have 
taken the costs of reusing urban residential or 
industrial water as a global constant of US$ 0.3 
per m?.9?? 


For each city-RCP combination, we calculated 
the minimum and maximum drought costs. We 
had to take a scenario approach here, as the 
literature contained no exclusions to a suitable or 
achievable mix of the three adaptation options. 
This probably varies by city, local situation 
and characteristics (for example, how much of 
the renewable surface water has already been 


consumed). Hence, for each city-RCP scenario, 
we calculated both the cost based on the 
cheapest solution and the cost based on the 
priciest solution to give a range for the potential 
cost of freshwater production. We also capped 
the reuse of urban residential/industrial water at 
50% of total drought volume, so it would only 
cover 50% of the total freshwater production 
needed to maintain business as usual. No global 
average or spatially varying number for this cap 
could be found in the literature, with examples 
ranging from 5% to 100% percent, depending 
on which residential/industrial water use was 
reused.** We thus chose a mid-way solution of 
50%. Estimates of desalination and reservoir 
capacity most often seemed capable of meeting 
the drought volumes in this study, so were not 
capped.» See Figures 5a and 5b for a flow chart 
summarising these steps. 
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Socio-economic analysis: 
Hydrological drought assumptions 


Even though drought volumes might be 
overestimated, drought costs are likely to be 
underestimated. We assume that the cost of water 
production/saving is a suitable proxy for actual 
drought cost and that the true drought costs are 
somewhere in between the minimum and maximum 
costs estimated here. Although this gives a 
reasonable estimate of the order of magnitude, the 
true costs are probably higher than the maximum 
costs estimated here due to chain effects within or 
outside the city’s economy. Actual drought costs - 
both direct and indirect - remain difficult to estimate, 
as the hazard is intangible. This makes it hard to 
discriminate between pure drought-related impacts 
and impacts that happen during but independently 
of the drought. 


Each water-production/-saving measure is assumed 
to be capable of fully solving each drought volume 
as estimated here, except for the reuse of urban 
residential/industrial water, which is capped at 50% 
of the total drought volume. Actual limits may differ 
from those assumed in this assessment due to high 
costs, physical production limitations or insufficient 
space for desalination facilities or dams, for instance. 
No further information on those limits was found and 
research on a local scale is needed to track down 
the corresponding numbers for this. Furthermore, 
desalination is assumed to be applicable anywhere, 
but could be relatively expensive for inland cities 
that have no saltwater bodies nearby, so need to 
transport water over large distances. 


Hydrological drought 

Yearly shortage of surface water 
in volume terms per 1,000 
citizens for each C40 city 


Unit costs of freshwater 
production/savings 
Unit costs of freshwater 
production/savings 
Unit cost estimates of freshwater 
production/savings using a mix 
of: 

A 


“a 


Increase in reservoir 
storage capacity 


Em 
$ 
a 
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Increase in 
desalination 
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é 4 


Reuse of urban 
residential/industrial water 


Figure 8: A visual summary of the methodology behind our assessment of socio-economic risk (cost to replace water deficit) from 


hydrological drought 
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Figure 9. Visual summary of the methodology behind our assessment of flood risk 


After calculating the socio-economic impacts for the 
various return periods, we estimated the expected 
annual impacts (EAIs), commonly referred to as the 
probabilistic risk. This step was only required for 
riverine and coastal flooding, as our hydrological 
drought calculations did not use return periods and 
already included an estimated annual damage value 
from the previous step. 


We calculated the EAI by integrating the impacts 
(the damages, exposed GDP or exposed population) 
over all of the different return periods, which 
effectively gave an estimate of the average impact 
per year (expected annual damage, expected 
number of people affected per year, expected 
effect on GDP per year). In the EAI, we included 
the impacts for each city for each possible GCM 


combination (riverine flooding only), RCP, SSP and 
impact type. In reality, there were many years with 
no impact and only one or a few years with a very 
high impact. The EAI indicators thus denote the 
average of all of these years. 


Standardising EAIs 

For each city, the EAI for riverine and coastal flooding 
and for hydrological drought is standardised per 
1,000 inhabitants. The reason for this is that total 
impacts are likely to be higher for larger cities, as 
these have greater populations, more buildings and 
a larger GDP (exposed elements) to be impacted. 
Part of this bias is taken away by standardising the 
EAIs for population. In other words, the standardised 
values denote the impact regardless of size. 


It is crucial to account for the fact that many regions of the world (mainly cities) are protected against flooding, for example, by 


dykes or levees. If these are ignored, as they are in most large-scale flood risk estimates, the risk will be grossly overestimated. 
To address this issue, IVM developed a database of global flood protection standards (FLOPROS),°° which we used in this study 
to exclude impacts for flood events below the current design standard for dykes and levees. It should be noted that FLOPROS 


includes three layers of protection standard based on three source categories: (1) the design layer, which holds information on an 


empirical basis and represents currently installed protection; (2) the policy layer, which fills gaps in the design layer where possible 


and holds information on policy regulations; (3) the model layer, which plugs any remaining gaps by modelling the protection 


standards in areas where no empirical or policy information is available.” Hence, FLOPROS includes protection standards that are 


already in place (design layer) and those that are desirable (policy and model layers). The latter may also be in place, but this could 


not be confirmed or denied from the data. Therefore, if there is a protection standard in place that guards the city against a flood 


event with a return period of 10 years, which equates to a 10% exceedance probability per year, all flood events with a return period 


of 1O years or less should not cause any damage to the city. 
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Definitions 


Climate risk 


CLIMATE 
RISK 


HARMFUL CLIMATE X 


R PROBABILITY OF A 
EVENT OR TREND 


Hazard x Exposure x Vulnerability 
The threatening x The elements x The resistance or 
= event (e.g., buildings, lack of resistance of 

(including its humans, the exposed 
probability and economic elements to the 
geographical extent) values) hazard 

present in 

affected areas 


Figure 10. A visual definition of climate risk 


and what is referred to as components of exposure 
and vulnerability. We follow the definitions of the 
United Nations Office for Disaster Risk Reduction 
(UNDRR) in this project. The following definitions 
of hazard, exposure, and vulnerability are from the 
UNDRR glossary,? along with a short addendum or 
explanation on each. 


The Intergovernmental Panel on Climate Change's 
(IPCC) Fifth Assessment Report (AR5) offers a 
conceptual framework for the assessment of risks 
associated with climate change? Climate risk is 
caused by harmful climate events that have negative 
impacts on cities worldwide. The consequences 
are the result of the interplay between the hazards 
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Hazard 

A process, phenomenon or human activity that 
may cause loss of life, injury or other health 
impacts, property damage, social and economic 
disruption or environmental degradation.9? It 
is a threatening event, including in terms of its 
probability and geographical extent. Note the 
word 'may', which implies that an event does 
not always cause (negative) impacts, as these 
depend on the exposure and vulnerability 
components. 


Climate risk 

The potential loss of life, injury or destroyed or 
damaged assets that could occur to a system, 
society or a community in a specific period of 
time, determined probabilistically as a function 
of (1) a climate-related hazard, (2) exposure, 
(3) vulnerability.55 Therefore, in facing a climate 
hazard, people, assets and natural resources 
may be exposed and vulnerable in different 
ways, depending, among other things, on their 
location and capacity to cope and adapt. The 
interplay of a hazard, the exposed elements and 
their vulnerability determines climate risk. 


General circulation models (GCMs): A model 
that predicts the climatic state of the earth based on a set 
of boundary conditions, namely solar radiation and GHG 
concentration. 


Shared socio-economic pathways (SSP): SSPs 
represent a set of scenarios - or plausible future 
worlds - that underpin climate change research 
and permit the integrated analysis of future 
climate impacts, vulnerabilities, adaptation and 
mitigation, including possible pathways for 
global society, economics and demographics. 
They can be categorised by the degree to which 
the different scenarios represent challenges to 
mitigation (reducing the sources or enhancing 
the sinks of GHGs) and societal adaptation to 
climate change. 


Representative concentrated pathway: RCPs 
aretrajectories of greenhouse gas concentrations 
resulting from human activity corresponding to 
a specific level of radiative forcing in 2100. 


Return period: For example, if a flood of volume 
X happens statistically once every 100 years, the 
return period of that flood is said to be 100 years. 


Rootzone moisture: The volume of moisture in 
the top 1 metre of soil in a given area. 


Exposure 

The situation of people, infrastructure, housing, 
production capacities and other tangible human 
assets located in hazard-prone areas.?' These 
are all the elements present in affected areas, 
such as citizens, flora and fauna, buildings and 
(critical) infrastructure. 


Vulnerability 

The conditions determined by physical, social, 
economic and environmental factors or processes 
which increase the susceptibility of an individual, 
a community, assets or systems to the impacts 
of hazards.*? This is the (lack of) resistance of the 
exposed elements to the hazard. 


Water province: A water province is a spatial 
unit of watersheds that belongs to exactly one 
country and one watershed area. If a water 
province is shared by multiple cities, we assume 
that the total water shortage is split equally 
among these cities. Where a city intersects 
with more than one water province, we add up 
the water deficits of all corresponding water 
provinces. Thus, each city is assigned their 
share of water shortages from adjacent water 
provinces. 
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Annex 1. 
Possibilities and 
limitations 

of global-scale 
models 


Flooding 


Riverine flooding 

Global riverine flood risk assessments generally 
include a full suite of hazard, exposure 
and vulnerability components.^ They have 
sporadically incorporated hazard and exposure 
dynamically. Dynamic vulnerability has only 
been included once, however? In all other 
cases, static hazard/exposure/vulnerability data 
were used. The spatial resolution of global-scale 
analyses in the most recent studies has varied 
from 30 arcseconds (“) to 2.5 arcminutes ¢‘).° 
Hazard has often been expressed as the depth 
and extent of inundation, population exposure, 
GDP or land use and vulnerability by depth- 
damage function.” To arrive at estimated flood 
areas and flood maps, most global models do 
not use a full hydrodynamic modelling scheme, 
as this requires high-resolution topographical 
features that are not yet available and have a 
high computational cost.9? Rather, like us, they 
compute discharge and use a volume-spreading 
approach to determine flooded areas. 


Coastal flooding 

Global coastal flood risk assessments are 
similar to riverine ones. They also consider all 
risk components, sporadically incorporating 
dynamic hazard and exposure. However, 
dynamic vulnerability has not yet been included. 
Spatial resolutions deviate slightly from riverine 
flood risk assessments too, varying between 
coastline segments and gridded fields of 30 
arcseconds (30”). Still, the risk components 
are explained using the same parameters as 
for riverine flooding. Although sea levels are 
calculated using a fully hydrodynamic model, 
the flood levels on land are not - a challenge 
similar to that discussed for riverine flooding.9? 
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Stormwater flooding 

Assessing stormwater flooding on a global 
scale is difficult, as it is highly dependent on 
local factors. To the best of our knowledge, 
the scientific literature does not contain any 
examples of global-scale stormwater flood risk 
assessments. It is most commonly assessed 
using flood models for a small area (such as city 
or even part of a city) that generate data on the 
depth and velocity of surface water associated 
with rainfall events of varying intensity. Guerreiro 
et al. (2017) have developed a continental-scale 
modelling approach to assess the stormwater 
flood hazard for 571 cities in Europe. The paper 
outlines some of the key challenges of such a 
continental approach, which would be amplified 
in the case of global-scale application. These 
include difficulties in obtaining the required 
hourly rainfall records, the low resolution of 
continental to global digital elevation models 
(DEMs) compared with those typically used for 
stormwater flood models, and a lack of data 
on local sewer systems, building shapes and 
infiltration of local green spaces. Opportunities 
to collect such data lie in high-resolution remote 
sensing and data science,” but also in local 
crowdsourced methods, such as community 
mapping," which is a promising avenue, 
particularly in strongly growing urban centres 
in developing countries. Guerreiro et al. (2017) 
demonstrate that current modelling capabilities 
and the requisite computing power make large- 
scale stormwater flood hazard assessment a 
possibility, provided these data challenges can 
be overcome. However, to include stormwater 
flooding in this research, we use the change in 
extreme rainfall as a proxy (see section 5.4.2). 


Due to their complex and diverse character, it 
is often difficult to find a comprehensive set of 
indicators for drought risk. This has led toa variety 
of approaches and datasets in global drought risk 
assessments. According to Ward et al. (2020), 
these assessments generally include indicators 
for the hazard itself and the corresponding 
exposure. Vulnerability has sporadically been 
included, but not yet dynamically. The spatial 
resolution of global drought risk assessments 
has been 0.2-0.5 degrees. As summarised by 
Ward et al. (2020), drought hazards have been 
expressed through a multitude of different 
indices, such as the Standardized Precipitation 
Index (SPD, Standardized Precipitation and 


Evapotranspiration Index (SPEI), Standardized 
Runoff Index (SRI), Palmer Drought Severity Index 
(PDSI) and the Weighted Anomaly of Standardized 
Precipitation (WASP). Exposure has been expressed 
as the effect on the population, GDP, road density, 
agricultural area and/or crop yields. Vulnerability 
is represented by the proportion of area irrigated 
per country or by proxies of social, economic and 
infrastructural vulnerability (for example, GDP per 
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capita, government effectiveness or percentage 
of irrigated agricultural land).’* Lastly, the risk is 
expressed using indicators such as maize yield 
losses or agricultural area or population affected.” 


Qualitative hydrological 
drought methodology 


The above methodology resulted in several C40 
cities” lacking a hydrological analysis for one of the 
following reasons: 


* According to the CWM database, the city 
only had one groundwater source. 

e The city was present in the database, but 
had no withdrawal point to link it to 
discharge. 

The city was not in the CWM database. 

The city was in the CWM database, but there 
was no information on its surface water 
sources. 


Consequently, a methodology had to be developed 
to qualitatively assess the impacts of hydrological 
drought on these cities. The cities for which data 
were missing are included in Table 2. 


To assess the relative impacts and benefits 
associated with the availability of urban water 
resources and their risk of hydrological drought, 
we used a multi-criteria analysis, comparing five 
indicators and evaluating different combinations. 


Five indicators were used: (1) water supply stress 
(WSS), (2) water-source type distribution (WST), 
(3) number of sources (S) (4) CDP drought 
intensity (CDP) and (5) recent drought activity (DA) 
(see details in Table 7). By assigning weights? to 
each indicator, it was possible to assess the drought 
impact in terms of these indicators. With each 
indicator carrying a maximum score of three, the 
total maximum score a city could have was 15.7 
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We conducted two multi-criteria analyses using 
these indicators. The first included all the indicators 
with the exception of cities that did not have any 
CDP reporting. 
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Both analyses had a minimum indicative percentage 
of 33% and a maximum percentage of 100%, due to 
the weighting of each indicator. For the final score, 
an average of the results of the two analyses was 
used for each city. 


To verify this methodology, it was applied to a 
regionally representative group of cities for which 
IVM had data.” As the normalised percentages 


Indicative 
percentage 


C40 ranking 
Stockholm 
Accra 
Nairobi 
New York 
Jakarta 
Beijing 
Delhi 
Sydney 


1 

2 
3 
4 
5 
6 
7 
8 


© 


London 


Rio de 
Janeiro 


Santiago 


y 
o 


Los Angeles 


Table 5. Ordinal hydrological drought ranking of cities 
The methodology was applied to all C40 cities for 
which data were missing. The results can be seen 


in Table 2. Using the indicative percentages for 
risk, we created a risk profile for all of the cities in 
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Population 


could not be quantitatively compared with the 
water-source volume difference for the 10-year 
return period, the ordinal ranking of the cities, from 
lowest risk to highest risk of drought, was used 
(see Table 5 ). Of the cities for which we had data, 
only three were not within 1 ordinal ranking when 
comparing the C40 and IVM rankings, namely, New 
York, Delhi and Santiago. 


IVM ranking 


Stockholm 


1,485,000 
3,013,000 
3,958,000 
8,468,000 
10,470,000 
18,079,000 
25,629,000 
4,844,000 


9,348,000 


12,380,000 
6,355,000 


14,081,000 


Accra 
Nairobi 
Jakarta 
Beijing 
Santiago 
Sydney 
London 


Rio de Janeiro 


New York City 
Delhi 


Los Angeles 


the analysis based on the Likert scale (see Table 
6). Cities were grouped into their corresponding 
categories based on this scale. 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 


Indicative percentage | Rank 
Zhenjiang 
Hangzhou 
Shenzhen 
Berlin 
Fuzhou 
Rome 
Madrid 
Dubai 
Kuala Lumpur 
Milan 
Singapore 
Medellín 
Stockholm 
Abidjan 
Quezon City 


Ekurhuleni 
Jaipur 
Miami 

Lima 
Accra 


Table 6. Cities’ indicative risk table based on the Likert Scale 


Water supply stress 


Water-source type 
distribution 


Number of sources 


CDP drought intensity 


Recent drought activity 


Indicator name 


Description 


Annual stress is defined as the ratio N/A 
of water withdrawn from a 
watershed to water available. By 
convention, a value above 0.4 is 
defined as stressed. The Urban 
Water Blueprint only evaluates 
surface water sources 


(missing 
information or 
primary sources 
are not surface 
water) 


Source mix, 
with no singular 
source type 
accounting for 
more than 50% 


The water-source typology 
categorised to include 
groundwater, surface water and 
desalination/other sources 


The total number of water sources City has more 
from which the city abstracts its than four water 
water sources 


A combined metric of probability Cities with a 


low drought 
intensity 


and magnitudetii of drought 
reported to the CDP by cities 


Drought[v] reported within a given 
city within the last 15 years 


No drought 
reported 


Table 7. Criteria used for the evaluation 
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Beijing 
Curitiba 
Delhi 
Jakarta 
Nairobi 
New York 
Shanghai 


Indicative percentage 


Washington, D.C. 
Johannesburg 


London 
Sydney 


Cape Town 
Melbourne 
Rio 


Copenhagen 


Hanoi 


Durban (eThekwini) 


Seoul 
Los Angeles 
Santiago 
New Orleans 


description 
Source 


(a value below 
or equal to 0.4) 


More than 50% 
of water 
sources are 
groundwater 


City has fewer 
than four water 
sources 


Cities with a 
medium drought 
intensity 


Single one-year 
drought 
reported 


(a value above 
0.4) 


More than 50% 
of water 
sources are 
surface water 


City has a 
single water 
source 


Cities with a 
high drought 
intensity 


Multiple or multi- 


year drought 
reported 


McDonald and Shemie 


(2014) 


CDP reporting by cities 


Multiple sources 
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Annex 3. 
Data sources 


Inter-Sectoral Impact Model 
Intercomparison Project 


(ISIMIP Fast Track) 


WATCH WFDEI 


City Water Map (CWM) 
Database 


Spatial Production 
Allocation Mode 
(MapSPAM) 


2UP 


Global Human Settlement - 
Settlement Model Grid (GHS- 
SMOD) 


Water provinces 


Natural Earth 10m coastlines 


Shuttle Radar Topography 
Mission Digital Terrain 
Model 90m V4.1 


Depth-damage curves 


FLOod PROtection 
Standards (FLOPROS) 
Database 


Multi-Error-Removed 
Improved-Terrain (MERIT) 
Digital Elevation Model 
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Description 


Bias-corrected global gridded climate data. 


Global gridded climate data for the 1901-2012 period, 
based on reanalysis datasets. 


Database that holds information on water extraction 
points and, wherever possible, withdrawal volumes for 
534 cities globally, based on web searches. 


Gridded map of the distribution of physical agricultural 
areas of 42 common crop types, based on a variety of 
input data. 


Global simulation of future urban growth, population 
distributions, and sub-national GDP developments, all 
based on SSPs. 


Global shapefile of city borders based on (1) population 
totals, (2) population densities and (3) build-up 
densities. 


Shapefile of intersections between administrative 
country borders and watershed borders. Each water 
province thus lies solely in one country and one 
watershed. 


Global vector dataset of coastlines based on World 
Data Bank 2 and NASA Mosaic of Antarctica. 


Relatively high-resolution global map of elevation 
based on satellite imaging. 


These curves depict the relationship between a certain 
inundation depth (water level) and the expected costs 
for acertain type of building and/or for a certain 


https://www.isimip.org/ 


https://www.isimip.org/ 


https://knb.ecoinformatics.org/view/doi963A10.5063962FF 1J 
67DWR 


https://www.mapspam.info/ 


https://www.pbl.nl/sites/default/files/downloads/pbl-2018- 
Towards-an-urban-preview_3255.pdf 


https://ghsl.jrc.ec.europa.eu/ghs_smod2019.php 


httos://www.sciencedirect.com/science/article/pii/S13 
64815218311447?casa_token=CgeXCjZv3NgAAAAA:u 
w697 pw4MGjn9- 

vr42bBMqkK9EZ Qw4Zeur3i_l¥dScDFvClfgdaZISrytE2h 
mJUJ5pZjfExILDJE 


httos://www.naturalearthdata.com/downloads/10m- 
physical-vectors/10m-coastline/ 


https://cgiarcsi.community/data/srtm-90m-digital- 
elevation-database-v4-1/ 


https://publications.jrc.ec.europa.eu/repository/handle 


location. These ones are construction cost surveys from /JRC105688 


multinational construction companies and are specified 
for each individual country. 


FLOPROS provides information on flood protection 
standards on a sub-national level (in this research). It is 
based on (1) empirical information where available, (2) 
policy regulations wherever available and where no 
empirical information is available, (3) modelling results 
for all other places. 


https://research.vu.nl/en/publications/flopros-an- 
evolving-global-database -of-flood-protection-standards 


http://hydro.iis.u-tokyo.ac.jp/-yamadai/MERIT. DEM/ 


PCR-GLOBWB 


(Geophysical Fluid 
Dynamics Laboratory - 
Earth Systems Model) 


GFDL-ESM2M 


HadGEM2-ES 


IPSL-CM5A-LR 


MIROC-ESM-CHEM 


NorESM1-M 
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A global hydrological model developed by Van Beek 
and Bierkens (2009) and Van Beek et al (2011). It is 
forced with gridded climate data. For the current period 
(1960-99), the model is forced with EU-WATCH forcing 
data (Weedon et al., 2011). For the future, it is forced 
using bias-corrected 


One of the GCMs used to capture the complex 
dynamics on Earth that make up the climate. Each GCM 
tries other sets of assumptions and methods to perform 
their calculations. As there is no ‘best’ model, scientists 
usually aim to use a couple of GCMs and take an 
average or median value from them to get a robust 
outcome of the future climate. 


One of the GCMs used to capture the complex 
dynamics on Earth that make up the climate. Each GCM 
tries other sets of assumptions and methods to perform 
their calculations. As there is no ‘best’ model, scientists 
usually aim to use a couple of GCMs and take an 
average or median value from them to get a robust 
outcome of the future climate. 


One of the GCMs used to capture the complex 
dynamics on Earth that make up the climate. Each GCM 
tries other sets of assumptions and methods to perform 
their calculations. As there is no ‘best’ model, scientists 
usually aim to use a couple of GCMs and take an 
average or median value from them to get a robust 
outcome of the future climate. 


One of the GCMs used to capture the complex 
dynamics on Earth that make up the climate. Each GCM 
tries other sets of assumptions and methods to perform 
their calculations. As there is no ‘best’ model, scientists 
usually aim to use a couple of GCMs and take an 
average or median value from them to get a robust 
outcome of the future climate. 


One of the GCMs used to capture the complex 
dynamics on Earth that make up the climate. Each GCM 
tries other sets of assumptions and methods to perform 
their calculations. As there is no ‘best’ model, scientists 
usually aim to use a couple of GCMs and take an 
average or median value from them to get a robust 
outcome of the future climate. 


https://github.com/UU-Hydro/PCR-GLOBWB_model 


https://www.gfdl.noaa.gov/earth-sy stem-model/ 


https://www.metoffice.gov.uk/research/approach/mod 
elling-systems/unified-model/climate- 
models/hadgem2 


https://cmc.ipsl.fr/international-projects/cmip5/ 


https://catalogue.ceda.ac.uk/uuid/d90ca0077e3344c 
7840ca56e49f89ee7 


https://gmd.copernicus.org/articles/6/687/2013/gmd- 
6-687-2013.html 
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Annex 4. 
Resolution tables 


Throughout our analysis we have had to used varying models 
with differing levels of resolution. We summarise the resolution 
of all datasets below. 


Climate hazard and 
socioeconomic analysis 


Indicator Resolution 


30” x 30” (around 1 km x 1 km at the 


Riverine floodin Flood volume per pixel 
9 PETE equator) 


The model resolution is 2.5 km from 


Coastal flooding Flood volume per pixel the coast (1.25 km in Europe). 


Change in the frequency of a one-in-10 year 


SPEM. 0.50 x 0.50 
precipitation event 


Stormwater flooding 


Population and GDP 30” x 30” 


Urban land use indicator from 
the 2UP model. The model has a 
fixed mix of three urban 
categories — residential, 2UP resolution = 30”x 30” 
commercial and industrial - and 
each class is assigned its own 
maximum potential damage. 
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